Spin-Hall effects in a Josephson contact 
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The Josephson tunnehng through a 2D normal contact with the spin-orbit split conduction band 
has been studied in the diffusive regime. Linearized Usadel equations for triplet components of 
the pairing function revealed a striking similarity to the equations of spin diffusion driven by the 
electric field in normal metals. Consequently, we predict that the out-of-plane spin-Ifall polarization 
accumulates towards lateral sample edges and the in-plane polarization is finite throughout the 
entire normal region. At the same time, the spin-Hall current is absent in the considered case of the 
stationary Josephson effect. 
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In connection with various spintronic applications, 
much interest have been attracted recently to spin-orbit 
interaction (SOI) effects on electron transport in normal 
metals and semiconductors. This interaction gives rise 
to fundamental transport phenomena, such as the spin- 
Hall effect (SHE) (for a review see jlj) and electric spin 
orientation [l|, 0] ■ These effects represent a direct man- 
ifestation of the spin-orbit coupling between spin and 
charge degrees of freedom in electron transport. At the 
same time, spin-orbit effects were also discussed for su- 
perconductors. Some works dealt with SFS junctions JSj 
(F stands for ferromagnet), others considered SNS [j], 
SN [6j] systems, or bulk superconductors [1, 0]. As was 
pointed out in Ref. 0, 0], SOI leads to admixture of 
triplet components to the pairing function. This sort of 
singlet-triplet coupling looks similar to the spin-charge 
coupling in normal systems. Therefore, one would ex- 
pect that phenomena closely related to SHE could mani- 
fest themselves in superconductors. At the first sight on 
this problem it becomes clear that, at least in the case 
of zero voltage across the junction, the spin-Hall current 
can not be generated as a linear response to the super- 
conducting current. The reason is that these currents 
have opposite parities with respect to the time inversion, 
while they must be equal in the stationary nondissipative 
superconducting transport. On the other hand, besides 
the spin currents, in normal systems SHE leads to spin 
accumulation near sample edges. Therefore, it is inter- 
esting to find out, if similar accumulation of magnetiza- 
tion takes place in superconducting systems. It should 
be noted that, despite formal similarities, such a magne- 
tization is fundamentally distinct from that induced by 
the normal SHE, since it is not subject to the energy dis- 
sipation accompanying spin diffusion and relaxation in 
normal systems. 

We will consider SHE and the electric spin orientation 
for a Josephson tunneling through a 2D normal contact 
(see Fig 1). The SOI there is represented by the Hamilto- 
nian Hso — cr • hfe, where cr is a vector consisting of Pauli 
matrices. The spin-orbit field hfc, which is a function 



of the electron wave vector fc, can be given, for exam- 
ple, by Rashba 0], or Dresselhaus [1| SOI, as well as by 
their combination. In this case the vector lies in the 
plane of the 2D system. The electron transport through 
the contact will be treated within the diffusion approx- 
imation, so that the length of the junction L, the elec- 
tron coherence length Lc and the spin precession length 
Lso — vp/h, where vp is the Fermi velocity and h is the 
angular averaged spin orbit field, are assumed to be much 
larger than the electron mean free path I. The electric 
voltage across the junction is set to zero. Hence, the su- 
percurrent is provided by the phase difference between 
two electrodes. The analysis of such a problem will be 
performed within a standard semiclassical treatment of 
Gor'kov's equations in the diffusion approximation (for 
a review see Q). Our goal is to derive linearized Usadel 
type equations and calculate the spin density induced by 
SHE. 

As far as the thermal equilibrium state is consid- 
ered, all observables of interest can be expressed via re- 
tarded and advanced Green functions. The correspond- 
ing Gor'kov's equations in the Nambu representation 
have the form 

(^^^ -^H- E''/'^^ G"'/''(A:, X') = 5{X - X') , (1) 
where r, a denote retarded or advanced functions, X = 




FIG. 1: Josephson contact with S and N denoting supercon- 
ducting and normal regions. 
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with the momentum operator k = —id/dr and the chem- 
ical potential /i. After averaging of initial Green func- 
tions over random positions of short-range impurities, the 
self-energy in (jT]) takes the form [lOj 
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G^/'^{t,t',r,r)TsS{r-r'), (3) 



where r is the elastic scattering time. Unperturbed 
Green functions are easily obtained from Eq.((T]). In 
the momentum representation and after the time Fourier 
transform they can be written as 

GO'-/"(a;, k) = {uj- t^E^ - <t • ± iVy^ , (4) 

where Ek = (fc^/2m*) — /i. Below we will perform calcu- 
lations for retarded functions and drop the labels r, a. 

Proximity to superconducting contacts results in an 
admixture to the Green functions of anomalous (propor- 
tional to Ti and T2) components. Also, these functions 
become inhomogeneous in space. In order to calculate 
them, we will follow a well known procedure in the frame- 
work of the semiclassical approximation fvi\. First, we 
perform the Fourier transform with respect to X ~ X' 
introducing, accordingly, the frequency and wave vector 
variables, to and k. The center of mass variables will be 
remained intact and denoted as r. Since the problem is 
stationary, the corresponding center of mass time vari- 
able is absent. Taking into account that variations of G 
in the scale of the Fermi wave-length are small, Eq. ([1]) 
should be expanded in terms of gradients d/dr. The next 
step is to simplify the self-energy part of Eq. ([1]) keeping 
there only terms linear in the anomalous part. Such a 
linearization can be done if the transparency of the SN 
contact is small, or the leads are close to the supercon- 
ducting critical temperature. Taking the sum of Eq. ([1]) 
and its conjugate one, and making use of the fact that 
hfc is an odd function of fc, for the anomalous part G12 
we obtain the equation 



{2Lo-v-ci+-)Gi2 - {hfc-cr,Gi2}- 

T 

2 [(^hfe^q • cr, G12] — Isc , 
where Shk q — (q • Vfc)hfe with q = —id/dr, and 
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(G?i.gi2 + gi2G22 



(5) 
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The lower labels in G" denote the matrix elements in 
the Nambu space and gi2 = J2k ^12- Usually, the quan- 
tum kinetic equation, such as (O, can be reduced to the 
Eilenberger l^] equation by integration with respect to 



the electron energy. In our case this procedure is not 
convenient because of electron energy spin splitting. In- 
stead, within the diffusion approximation, from Eq. ^ 
we will express G12 in terms of 312, and taking its sum 
over k obtain the closed diffusion equation for gi2. Be- 
fore doing this, we transform the 2x2 matrix Gi2,q/3 to 
the conventional pairing function F^-^ = Gi2,a/3, where /? 
denotes the spin projection opposite to /?. Further, it is 
convenient to decompose F into triplet Fi,F_i,Fo and 
singlet Fs components as 



_ F12 + F21 _ F12 — F21 

Fi = Fii , F_i — F22. 



(7) 



The corresponding density function / = F will also 
be represented in a similar way. After this transforma- 
tion, it is easy to see that the last term in the l.h.s. of 
([5]) is responsible for a coupling between the singlet and 
triplet components of the pairing function. Besides, the 
singlet-triplet coupling also originates from the spin de- 
pendent parts of Gil ^^'^ ^22 Eq. ([6]). Due to such 
coupling, the triplet component of F is generated within 
the junction between two singlet superconductors. 

For simplicity, when deriving the diffusion equation, let 
us assume that SOI is strong enough, so that Lso <SC Lc- 
Further, considering I^c together with the last term in 
the l.h.s. of ^ as sources, we resolve Eq.® performing 
expansion in [v ■ q)T and h^r up to the second order. 
Finally, we obtain the following diffusion equation for 
the triplet pairing function /„j — (i/nNp) Fm, {m = 
0,1,-1): 
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dr 



2J.hfc >/ + M/, 



(8) 



where J is the vector of 3 x 3 angular moment operators 
and (...) denotes the angular averaging over the Fermi 
surface. The triplet-singlet coupling is given by 

Mo = , M±i = ^ (hj (hfc X 6hkA))fs , (9) 



with hj 



The singlet fs satisfies the usual 
Usadel equation Q with an additional term which is Her- 
mitian conjugate to M. Since this term is small, we will 
neglect a corresponding correction to fs in Eq.([8]). Hence, 
fs is given by the well known unperturbed solution in the 
SNS contact. Since it varies within the scale Lc ^ iso, 
we neglected all contributions to M with higher powers 
of gradients, as well as terms proportional to lu ^ D / L^, 
where D is the diffusion constant. 

Without the last term in the r.h.s., Eq.® formally 
coincides with the spin diffusion equation for 2DEG in a 
zero electric field [13| . The spin diffusion equation in the 
presence of the electric field has been derived in Ref. [l| 
for the case of the Rashba SOI, and for a general SOI in 
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|15j . After a linear transformation to spin density 
variables Eq.® will also coincide with these equations, 
if, apart from a constant factor, /g is formally identified 
with the electric field potential. Hence, a coupling of the 
spin to the electric field in normal spin transport appears 
to be very similar to the singlet-triplet coupling in Eq.®. 

Let us consider an example of the Rashba SOI. In this 
case = aky and /i| = —akx- For a homogeneous in 
y-direction case all functions depend only on x and we 
get /o = 0, /i = with /i satisfying the equation 
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is the D'yakonov-Perel' spin re- 



where Fso = 2Ta'^k'j 

laxation time [3]. The small l.h.s. of Eq.® has been 
neglected in ([TO)) . Boundary conditions at x — ±L/2 
can be written in a way similar to a singlet SN interface 
17| . At least in the linearized approximation the bound- 
ary conditions contain only characteristics of one-particle 
transmission. Therefore, they can be easily generalized 
to the case of a triplet pairing. Following calculations of 
Ref. [17] we obtain 



(/is - ff/ijv) 



x=^L/2 



= ±b 



Of: 
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(11) 



where the labels S and N denote superconductor and 
normal sides of SN contacts at a: = ±L/2, and g — 
\uj\/y^{uj + i0+)2 — |Ap is a DOS factor for a super- 
conductor. The characteristic length b depends on the 
SN barrier transmittance. For our choice of parameters 
b 3> Lso- The same equation pT]) takes place for fg. 
At the low SN barrier transmission one may use the so 
called rigid boundary conditions and set /is — 0. At 
the same time, the singlet paring function /ssI=fL/2 — 
gAexTp{zLi(j>)/uj. Neglecting the third derivative of fg, 
the solution of Eq. PH)) can be written as 



(12) 



where ip{x) is a linear combination of exp(±fca;), with 
k = y/D/r = 1/ Lso- It is easy to see from (fTTj) that at 
A:6 ^ 1 the first term dominates in Eq. (fT2)) . Therefore, 
ijj will be neglected below. 

Our next step is to calculate the spin polarization den- 
sity associated with triplet components of the pairing 
function. This polarization is given by 



duj 
2^ 



np{Lo) X 



Tr[a\Gl,,{Lo,r)~Gl,^{uj,r)], (13) 

where np is the equilibrium Fermi distribution function. 
It is easy to see that the nonzero value of Eq. (fT3|) 
is provided by triplet components of anomalous Green 
functions which contribute to Gn with a correction term 



cx p. Up to the leading second order with respect to 
fs and keeping only the linear terms of the triplet fm 
[m = 1,-1,0), for the retarded function we obtain from 
Eqs. (fTTO 



r/ai _ Tl fiS'' 



ttNf 2 



l_^fr/af+r/a^jr/aj+r/a^^^ (14) 



where fy = (/i + /_i)/2 and = -z(/i - /_i)/2. The 
conjugate functions f^{uj) = — /*(— a;). 

In the case of Rashba SOI fx ^ fo = and fy — fi. 
The latter is given by Eq. (HI]). Then, from (HI]) it 
immediately follows that only the y-projection of the spin 
density is finite. Using the relations /"(w) = /J(— tj) and 
f^{uj) = —f^j^{—Lj) (m = 1,-1,0), wc arrive to the spin 
polarization 



Sy{x) = eNpar 



Jix) 



(15) 



where adc is the dc conductivity of the normal metal and 
J is the Josephson current density 



J = 



eD 



pj d^npi.m:J^-J^fr)-ir^a)]. 

(16) 



The spin polarization (jlSp coincides with polarization in- 
duced in normal metals by the electric field Q, if the 
Josephson current is substituted for the normal dissipa- 
tive dc current Jdc = <^dcE. It is easy to check that this 
analogy takes place also for the Dresselhaus SOL with 
a little more complicated expression for 5'*(x) [15|. An 
important distinction from the electric spin orientation 
in normal metals is that due to the charge neutrality, 
Jdc — const in the x direction, while the supercurrent 
varies inside the contact. Similar effect has been pre- 
dicted by Edelstein Q for bulk superconductors and at 
NS boundary, providing the supercurrent flows along the 
SN interface. 

Let us now check, if the analogy with the electric spin 
orientation extends to the spin-Hall effect. Hence, our 
goal is to calculate Jy, which is the y projection of a 
spin flux polarized in the z-direction. The corresponding 
spin current operator can be written as 3y = {az,^y\ 
where the velocity Vy = ky/m* -\- d{cr ■ hk)/dky. Since it 
has been assumed that hz = one gets = azky/m*. 
The spin-Hall current JsH, in its turn, can be derived 
from Eq. p^ . with cr' substituted for J^. Keeping the 
same leading terms as in calculation of the spin density, 
we arrive to Jsh — 0. This result does not depend on 
whether hp, is given by the Rashba or Dresselhaus inter- 
actions. That is very distinct from the normal spin-Hall 
effect, where in the diffusive regime the spin-Hall conduc- 
tance is zero for the Rashba SOI, but finite for the cubic 
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Dresselhaus interaction [l8|. In general, as it was dis- 
cussed above, the zero value of JsH in superconducting 
transport follows from the time inversion symmetry. 

Besides JsH, in normal systems the DC current to- 
gether with SOI gives rise to accumulation of the z- 
component of spin at the lateral edges of the sample 
[H in, HO]. In the case of the Josephson junction the 
0-projection of the spin density is given by Eq. (fT3ll and 
the first term in Eq. (fT4|) . Hence, it is proportional to the 
/o component of the pairing function which, in its turn, 
can be found from Eq. ([5]) . For simplicity, let us consider 
hard wall boundaries of 2DEG dX y = ±Ly/2. In this 
case one can borrow the boundary conditions for Eq. ([8|) 
from Ref. 15, [ISl- In normal systems these conditions 



correspond to the vanishing spin current at y = ±Ly/2. [4 
In our case similar equations can be written for triplet 
"currents" j = "^^^vF. We thus have j^|y=±L /2 = 0, 
where the triplet component is given by 

= -D^ - 2tT{{vy[hk X {V2f + rShk^eifs)]) ■ (17) 

The first term in this equation is the diffusive current, 
the first term in the brackets is determined by the spin 
precession in the effective spin-orbit field, and the last 
term looks as the spin-Hall current in the normal spin 
transport. As far as fg is treated a slowly varying func- 
tion of X, thus allowing one to ignore its higher gradi- 
ents together with edge terms like tp in (fT2|l . the analysis 
of Eq. ^ with the above boundary conditions is the 
same, as for SHE in normal systems. Henceforth, follow- 
ing Ref.0,12^ one may conclude that /o = for Rashba 
SOI, but fo is finite in the case of the cubic Dresselhaus 
interaction. From Eqs. p^ and (fT4|l it is immediately 
seen that in the former case ~ 0. For the Dressel- 
haus SOI the solution of Eq. ([5]) has the form fo = 
x{y){dfs/dx), where x is a real odd function of y. Then, 
Eqs. 111,(1141) and CID give = -eNFx{v){J / (^dc)- The 
fimction in its turn, have been calculated in Ref. [l5| . 

In conclusion, the spin-Hall effect induced by a super- 
current across an SNS junction has been studied in the 
diffusive regime for a relatively strong {Lso 3> Lc) SOI in 
the 2D junction and for low conducting SN barriers. We 
found out that, although the spin-Hall current is forbid- 
den by the time inversion symmetry, in the case of cubic 
Dresselhaus SOI the out-of-plane magnetization is accu- 
mulated near sample edges at y = ±Ly/2, in a very close 
analogy to SHE in normal systems. Also, similar to the 
electric spin orientation, the spin polarization parallel to 
2DEG is finite throughout the entire N-region. 
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